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Abstract. We study AAK-type meromorphic approximants to functions of the 
form 



where R is a rational function and A is a complex measure with compact reg- 
ular support included in (—1,1), whose argument has bounded variation on 
the support. The approximation is understood in L p -norm of the unit circle, 
p > 2. We dwell on the fact that the denominators of such approximants sat- 
isfy certain non-Hermitian orthogonal relations with varying weights. They 
resemble the orthogonality relations that arise in the study of multipoint Pade 
approximants. However, the varying part of the weight implicitly depends on 
the orthogonal polynomials themselves, which constitutes the main novelty 
and the main difficulty of the undertaken analysis. We obtain that the count- 
ing measures of poles of the approximants converge to the Green equilibrium 
distribution on the support of A relative to the unit disk, that the approxi- 
mants themselves converge in capacity to F, and that the poles of R attract 
at least as many poles of the approximants as their multiplicity and not much 
more. 
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1. Introduction 

This paper is concerned with the asymptotic behavior of certain meromorphic 
approximants to functions of the form 





(1.1) 
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where R is a rational function, holomorphic at infinity, and A is a complex measure 
compactly and regularly supported on (—1,1). 

The meromorphic approximants that we consider are optimal, for fixed num- 
ber of poles in the unit disk, with respect to an L p -norm on the unit circle. When 
studying them, we assume that supp(A) and all poles of R lie in the open unit 
disk, so that F is indeed p-summable on the unit circle. The asymptotics are then 
understood when the number of poles grows large. In the case where p = oo, this 
type of approximant was introduced by V. M. Adamyan, D. Z. Arov, and M. G. 
Krcin in their famous paper [T]. Here we deal with their natural generalization 
to L p , although we restrict ourselves to the range 2 < p < oo for technical rea- 
sons to be explained later. The meromorphic approximation problem also has a 
conformally invariant formulation on Jordan domains with rectifiablc boundary, 
to which the results of the present paper transpose with obvious modifications if 
supp(A) is contained in a closed hyperbolic geodesic arc rather than a segment. 
The interested reader will have no difficulty to carry out this generalization using 
the construction of [7] Sec. 5]. 

The study of best meromorphic approximants is quite recent. After the de- 
velopment of the Adamyan- Arov-Krein theory for p = oo in pQ, the latter was 
extended to the range 1 < p < oo by F. Seyfert and the first author in |10j , and 
independently by V. A. Prokhorov in [30]. However, it is only for p > 2 that the 
authors of [10] were able to express the error in terms of (generalized) singular- 
vectors of a Hankel operator and subsequently to obtain integral formulas for that 
error when the approximated function is represented as a Cauchy integral. These 
formulas make connection with non-Hermitian orthogonality, and form the basis 
of the present approach. 

The AAK theory had considerable impact in rational approximation, for on 
retaining only the principal part of a best meromorphic approximant to a function 
analytic outside the disk and sufficiently smooth on the circle, one obtains a near- 
best rational approximant to that function |15j . This is instrumental in Parfenov's 
solution, for simply connected domains, to the Gonchar conjecture |17] on the 
degree of rational approximation to holomorphic functions on compact subsets of 
their domain of analyticitjj^] and also for instance in Poller's converse theorems on 
smoothness of functions from their error rates in rational approximation |25l I26| . 
The same principle is also at work in the articles H] that deal with rational 
approximation to Markov functions, that is, functions of the form (1.1 1 where 
R = and A is positive. 

Another connection between meromorphic approximation and rational ap- 
proximation that ought to be mentioned occurs when p — 2. In this best 
meromorphic approximant to a function analytic outside the unit disk is in fact 
rational, and turns out to be a special type of multipoint Pade approximant that 
interpolates the function with order 2 at the reflections of its poles across the 



1 the proof of this conjecture was later carried over to the multiply connected case by Prokhorov 
in [29]. 
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unit circle |23l |9] . Of course the interpolation points are not known a priori which 
accounts for the nonlinearity of the problem. Nevertheless, after the work in [15] . 
this connexion was used in [TT] to establish the convergence rate of best L 2 ratio- 
nal approximants to Markov functions. The forthcoming results will, in particular, 
generalize these results to a larger class of functions. 

Best meromorphic approximants to Markov functions were studied per se by 
E. B. Saff, V. Prokhorov and the first author in [5J. Using results from [1] to make 
connection with orthogonality, these authors prove (and give error rates for) the 
uniform convergence of such approximants, locally uniformly on C \ /, whenever 
1 < P < oo provided that A satisfies the Szego condition: log dX/dt e L 1 ^). 

The present paper appears to be the first to deal with convergence of best 
meromorphic approximants to general functions of the form (1.1) in the case where 
A is a complex measure (the special case when F has two branchpoints and no 
poles is treated in [TO] Sec. 10]). 

This paper is organized as follows. Section [2] deals with meromorphic ap- 
proximants that are solutions (more generally: critical points) of the meromorphic 
approximation problem for functions of the form (1.1). We discuss the asymp- 
totics of poles as being the limit zero distribution of polynomials satisfying certain 
non-Hermitian orthogonality relations with respect to varying measures. We ap- 
ply the results to the convergence in capacity of these approximants, and to the 
convergence of some of their poles to the polar singularities of F. All the proofs 
are presented in Section [3] Some computational results are adduced in Section [4] 
and the Appendix contains necessary material and notation from potential theory 
that we use throughout the paper. 

Finally, we mention that all the results below have their counterpart in diag- 
onal multipoint Pade approximation, where they allow more irregular A than can 
usually be handled via classical results (see [19]). We do not include this to keep 
the size of the paper within reasonable bounds, and refer the interested reader to 
[31 Ch. Ill] or P]. 



2. Meromorphic Approximation 

Let A be a complex Borel measure whose support S := supp(A) C (—1, 1) consists 
of infinitely many points. Denote by |A| the total variation measure. Clearly A is 
absolutely continuous with respect to |A|, and we shall assume that its Radon- 
Nikodym derivative (which is of unit modulus |A|-a.e.) is of bounded variation. In 
other words, A is of the form 



d\(t) = e lipW d\X\(t), 



(2.1) 
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for some real-valued argument function ip such tha' 



2 



N 

V(<p,S):=sup{Y,\<p(x j )-<p(x j - 1 )\} <oo, (2.2) 

where the supremum is taken over all finite sequences xq < x\ < . . . < xjy in S as 
N ranges over N. 

For convenience, we extend the definition of ip to the whole of R as follows. 
Let I := [a, b] be the convex hull of S. It is easy to see that if we interpolate ip 
linearly in each component of I\ S and if we set ip(x) := lim t ^ aj t es fit) f° r x < a 
and <p(x) := lim i _ 1 .(, ) t es <p{t) for x > b (the limits exist by (2.2)), the variation of <p 
will remain the same. In other words, we may arrange things so that the extension 
of ip, still denoted by ip, satisfies 

V(<p,S) = V(p,R) =: V(p). 



Among all complex Borel measures of type (2.1 )-(2.2 1, we shall consider only 
a subclass BVT defined as follows. We say that a complex measure A, supported 
on (—1,1), belongs to the class BVT if 

(1) supp(A) is a regular set; 

(2) there exist positive constants c and L such that, for any x £ supp(A) and 
S £ (0, 1), the total variation of /i satisfies |A|([a; — 5,x + 8]) > c5 L ; 

(3) A has an argument of bounded variation. 

Denote by V n the space of algebraic polynomials of degree at most n and by 
M. n the subset consisting of monic polynomials of degree n whose zeros lie in the 
open unit disk, D. 

Define 

F(z):= J d ^ + R s (z), (2.3) 
with A £ BVT and R s £ 1Z s -i. s , where 

fcm,n ■= {Pmfan ■ Pm G Vm, <?n G M n } 

is the set of rational functions of type (m, n) with all their poles in D. Hereafter 
we shall denote by Q s the denominator of R s , assumed to be in irreducible form, 
which is a monic polynomial with zeros in D of the form 

Qs(z)=H(z-r 1 r^, (2.4) 

where S' is the set of poles of R s and m(rf) stands for the multiplicity of rj £ S'. 
Thus, F is a meromorphic function in C \ S with poles at each point of S' and 
therefore holomorphic in C \ S, where 

S:= SUS'. 



2 Note that e llfi has bounded variation if and only if ip can be chosen of bounded variation. 
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Note in passing that F does not reduce to a rational function since S consists of 
infinitely many points, cf. |7j Sec. 5.1] for a detailed argument. 

In this paper we consider the behavior of certain meromorphic (AAK-type) 



approximants to a function F of the form (2.3 1. Of particular importance to us 
will be the asymptotic behavior of the poles of the above-mentioned approximants. 
The latter will be quantified in terms of the weak* convergence, when the number 
of poles increases indefinitely, of the counting measures of these poles. By defini- 
tion, the counting measure of the poles of a meromorphic function is the discrete 
probability measure with equal mass at each finite pole, counting multiplicities. 
The weak* convergence is understood in the usual sense where measures, endowed 
with the norm of total variation, are regarded as the dual space to continuous 
functions with compact support. 

We denote by H p , p £ [l,oo], the Hardy space of the unit disk consisting of 
holomorphic functions / such that 

||/H5:= sup i- / \f(rO\ p m < oo if j>€[l,oo), 

(2.5) 



0<r<l 27T Jf 



H/lloo := sup |/(z)| < oo if p = co. 

It is known ( |141 Thm. 1.5.3]) that a function in H p is uniquely determined by 
its trace (nontangential limit) on the unit circle, T, and that the L p -norm of 
this trace is equal to the H p -norm of the function, where LP is the space of p- 
summable functions on T. This way H p can be regarded as a closed subspace of 
L p . Analogously, we define Hq , p £ [1, oo ], co nsisting of holomorphic functions in 



' \ D that vanish at infinity and satisfy (2.5 1 this time with 1 < r < oo 



Now, the meromorphic approximants that we deal with are defined as follows. 
For p £ [1, oo] and n £ N, the class of meromorphic functions of degree n in L p is 

H p :=H p + n n _ 1>n = H p B-\ (2.6) 

which is a closed subset of LP (it is in fact weakly closed if 1 < p < oo and weak* 



closed ifp = oo, see [TO] Lemma 5.1.]). In (2.6) we denote by B n the set of Blaschke 



products of degree at most n, consisting of rational functions of the form 

,q(z) 



biz) = e lc , q£Mk, where q(z):=zq(l/z), k < n. 
q{z) 

We shall call q the reciprocal polynomial of q'mVk- We also say that b is normalized 
if e IC = 1. Thus, the members of B n are rational functions of degree at most n 
holomorphic in D and having modulus 1 everywhere on T. 

Our best-LP meromorphic approximation problem can now be stated as fol- 
lows. 

MA(p) : Given p £ [1, oo], / £ LP , and neN, find g n £ H p such that 

||/-0n||p= infj|/-. 9 || p . (2.7) 
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Originally this problem was solved for the case p = oo by V. M. Adamyan, 
D. Z. Arov, and M. G. Krein in [1] and the solution came through operator theory. 
The most accessible reference to this result is perhaps [35]. Later F. Seyfert and 
the first author [TO], and independently V. A. Prokhorov j3Uji generalized it to the 
case 1 < p < oo, but it is only in |10j and when p € [2, oo] that concrete equations 
were obtained for the approximants. These form the basis of our approach, and 
presently limit its scope to p > 2. 

This solution of MA(oo) is known to be unique, provided that / belongs 
to the Douglas algebra H°° + C(T), where C(T) denotes the space of continuous 
functions on T [T]. In particular, the solution to MA(oo) is unique when / is of 



type (2.3 1 since the latter is analytic in some neighborhood of the unit circle (i.e. 
in the complement of S). When p < oo, a solution needs not be unique even if / 
is very smooth |10l Sec. 5]. Therefore, when making a statement about a sequence 
{g n } of solutions to MA(p), it is understood that a particular solution has been 
selected for each n and that the statement holds true regardless the selection. 



Given F as in (2.3 1, we shall be interested in three types of questions: 

(a) What is the asymptotic distribution of the poles of best-L p meromorphic ap- 
proximants to F as n tends to oo? 

(b) Do some of these poles converge to the polar singularities of F? 

(c) What can be said about the convergence of such approximants to F? 

As it is the case of interest here, we shall restrict our discussion to the situa- 



tion where the approximated function is of the form (2.3 1, and accordingly we write 
F instead of /. We should note that MA(2) reduces to rational approximation. 
Indeed, L 2 can be decomposed into the orthogonal sum of H 2 and its orthogonal 
complement Hq , which consists of analytic functions in C \ D vanishing at oo with 
norm sup r>1 ||/(V-)I|2 (compare (2.5 1). Now, since F is the Cauchy transform of a 



measure supported in D, it belongs to Hq. Thus for any g = (h + p n -i/q n ) & H 2 
with h £ H 2 and p n -i/<ln & rational function in Hq, we get by orthogonality 

\\F- g \\l = \\h\\ 2 2 + \\F~p n ^/q n \\l 

Clearly then, for g to be a best approximant h must be zero. It also turns out in 
this case that best approximants interpolate F with order 2 at the reflections of 
their poles across T and also at infinity with order 1. Thus, one can regard MA (2) 
as an interpolation problem of the multipoint Pade type where the interpolation 
points are implicitly defined by the solution. Despite this, we shall not distinguish 
p = 2 from the other cases but rather keep a unified operator approach. 

Let us denote by P + and P^ the analytic and anti-analytic projections acting 
on Fourier series as 

(+oo \ +oo / +oo \ —1 

k— — oo / k—0 \k— — oo / k— — oo 

By a well-known theorem of M. Riesz [14], P+ : D> -► H p and P_ : W -► Hq 
are bounded when 1 < p < oo, and when p = 2 they are just the orthogonal 
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projections associated to the orthogonal decomposition: L 2 = H 2 Hq. When 
/ £ L 1 , we simply regard P+(/) and P_(/) as Fourier series of distributions. 
For each p G [2, oo], the Hankel operator with symbol F S V is given by 

Ap : Hp' -> H% 

u i — ► P_(Fu), 

where p' is conjugate to p modulo 2, i.e. l/p+ = 1/2. 

For n = 0, 1, 2, ... , the n-th singular number of the operator Ap is defined to 

be 

a n {Ap) := inf ||||Af — r|||, T : _ff p — > Hq a linear operator of rank < nj , 

where 1 1 1 • 1 1 1 stands for the operator norm; when p = 2 we assume in addition that 
r is weak* continuous. 

Note that if p = oo then p' = 2, hence Ap operates between Hilbert spaces, 
and since it is compact]^] the o- n (Ap) are just the singular values of Ap, that 
is, the square-roots of the eigenvalues of A* F Ap arranged in nonincreasing order; 
throughout A* F indicates the adjoint of Ap. When 2 < p < oo, the usual eigen- 
vector equation gets replaced by a nonlinear equation of Hammerstein type. More 
precisely, to each a n (Ap) there exists (at least one) v n £ H p of unit norm, whose 
inner factoi^jis a Blaschke product of degree at most n, such that 

A* F Ap(v n ) = a 2 n (A F )P + (\v n \P'- 2 v n ) if p > 2, 

A* F A F (v n ) =P+ (\A F (v n )\ 2 v n ) and \\A F (v n )\\2 = (Tn(A F ) if p = 2. 

Such a v n will be called a n-th singular vector for Ap. From the definition it follows 
that a n-th singular vector can be factored as 

V n = b n W n (2.9) 

where b n £ B n and w„ is an outer function of unit norm in HP . Actually, upon 



pairing with v n , equation (2.8) implies that v n £ B n if p — 2, so that w n = 1 in 
this case. 

The solutions to MA(p) turn out to be exactly the functions of the form 
(PES Thm. 8.2]) 

A F (v n ) P+(Fv n ) 
9n = F = , (2.10) 

V n V n 

where v n is some n-th singular vector for Ap. Moreover, we have that \\F — g n \\ p = 
a n (Ap). Such a g n is called a best meromorphic L p -approximant of order n to F. 



3 This can be deduced from the fact that F e C(T), see [23 Thm. 1.5.5]. 

4 Recall that any function h from a Hardy space can be written as the product of an inner 
and an outer factors. The outer factor of h is equal to the exponential of the Riesz-Herglotz 
transform of log \h\, and the inner factor, which is a H°° function unimodular a.e. on T, consists 
of a Blaschke product (finite or infinite) and a singular inner factor which is the exponential of 
the Riesz-Herglotz transform of a singular measure on the unit circle. More on the inner-outer 
factorization of H p functions can be found, for example, in 1141 . 
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The notion of a best approximant can be further weakened to the notion of a 
critical point. By definition, a function g n is a critical point of order n in MA(p) 
if and only if it assumes the form 

Ap(v n ) 

9n = F , (2.11) 

v n 

where v n is a H p function of unit norm (a Blaschkc product if p — 2) whose inner 
factor lies in B n \ B n -\ and which is such that 

A* F A F {v n ) = 7 »P+ (kl p '- 2 «„) , 7™eM, if p>2, 
A F A F (v n )=P+(\A F (v n )\ 2 v n ) if p = 2. 



The difference with (2.8 1 is that here j n (or ||^4i?(u„)||2 if p = 2) needs not be 
equal to a n (A F ). With a slight abuse of language, we will continue to say that 
v n is a singular vector associated to g n although 7„ may no longer be a singular 
value. Note that, as for best meromorphic approximants, v n reduces to a Blaschke 



product when p — 2. Thus v n has an inner-outer factorization of type (2.9 1 where 
b n has exact degree n and w n = 1 if p = 2. 

Although their definition is a little technical, critical points are just those 
g n G for which the derivative of \\F — g n \\p with respect to b n € B n and 



h £ i? p in factorization (2.6 1 does vanish [10] . Beyond best approximants, the most 
important critical points are local best approximant [101 Prop. 9.3.]. By definition, 
a local best approximant is some g n S H? for which there exists S > such that 

g e HP and ||.g - .g„|| p < <5 imply \\F - g n \\ p < \\F - g\\ p . 

The reason why we introduce critical points is that all a numerical search can yield 
in general is a local best approximant, and we feel it is important that our results 
should apply to computable objects. 

When p € [2,oo), best and local best approximants have exactly n poles, 
counting multiplicities, hence they are a fortiori critical points of order n ( |10[ 
Prop. 9.2 and 9.3]); such critical points are called irreducible. For p = oo the 
critical points are just the best meromorphic approximants, which are unique 
for each fixed n. So the notion is nothing new, but it may happen that a best 
approximant out of H^ has less than n poles. However, each time the number 
of poles of g n increases with n, it jumps to the maximum value n, in particular, 
there exists a subsequence of natural numbers, say No = No(-F), such that for each 
n E No the best approximant g n has exactly n poles in D, i.e. it is irreducible ( |10l 
p. 114]). Since the behavior of the poles of best approximants from is entirely 
characterized by this subsequence, hereafter we say "a sequence of irreducible 
critical points of order rt" to mean if p — oo that we pass to a subsequence if 
needed. 

The three theorems stated below constitute the main results of the paper. For 
the definitions of capacity, Green equilibrium distribution, and condenser (Green) 
capacity, the reader may want to consult the appendix. 
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Theorem 2.1. Let p G [2,oo], p' the conjugate exponent modulo 2, and {g n } n <£N be 
a sequence of irreducible critical points of order n of MA(p) for F, where F is 



given by [2.3)-{2.J i ) with X G BVT. Then the counting measures of the poles of 
g n converge to /i(s,T); the Green equilibrium distribution on S relative to D, in the 
weak* sense. 

The previous theorem gives one answer to question (a). The next one ad- 
dresses question (c) by stating that approximants behave rather nicely toward the 
approximated function, namely they converge in capacit^to FonD\ S, and in 
the case p = 2 uniformly in C \ D. Moreover, n-th root estimates for the error are 
provided. 



Theorem 2.2. Let F and {g n } n< £n be as in Theorem 2.1 Then 



K F- 9 „)wr»^xp{ar"M-^} (2.i3) 

on compact subsets ofH>\S, where U^ s ' r) is the Green potential of ^(s,t) relative 
to D and C —> denotes convergence in capacity. Ln addition, in the case of rational 
and AAK approximation (i.e. when p = 2 and p = oo, respectively) it holds that 

| ( ,_ fc)W |.^5„{-_^- 0l «-( 1 /»)} (2-14) 
on closed subsets o/C\(DUS*). Moreover, for p = 2, it also holds that 



limsup |(F - g n ) (z) 1 1 ' 2n < exp {- 

n— >oo { 



< (s ' T '(l/z-U (2.15) 



cap(5*, T) 

uniformly in C \ D, where S* is the reflection across T of S. 

As a consequence of the previous theorem, we can prove a result on best 



rational approximation function of the type (2.3 1 in C\D that is classical in scope 
(cf. |19l Thm. 1']), but new if A vanishes on a subset of positive capacity of the 
convex hull of S. In what follows || • \ \k stands for the supremum norm on a set K. 



Corollary 2.3. Let F be given by (2.3\)-l2.4\) with X G BVT and E := C \ D. Then 

1 



Km Pn (F,E) 1/2n = expj - 



cap(supp(A), T) 



(2.16) 



wh 



ere 



Pn (F,E):= inf ||F-r|| E . 

To approach question (6), we need to introduce some more notation. For any 
£ ^ G C, we let Arg(£) G (— tt,it] be the principal branch of the argument and 
for £ = we set Arg(O) = ir. With this definition, Arg(-) becomes a left continuous 



3 See the appendix for the definition of convergence in capacity. 
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function on R. Now, for any interval [a, b] C K we can define the angle in which 
this interval is seen at £ G C by 

Angled [a, b}) |Arg(a - - Arg(6 - 0|. 

We define additi vely this angle for a system of disj oint closed intervals : if { [a,j , b 3 ■] } ™ = 1 
is such a system, then the angle in which it is seen at £ is defined by 

m 

9(0 := Angled [a^bj]). (2.17) 
j=i 

Note that < 9(0 < ir and 8(0 = w if and only if £ G U[aj,6 J -]. The notation 
#(£) does not reflect the dependency on the system of intervals, but the latter will 
always be made clear. Further, for any point z G C define the lower and upper 
characteristics m(z),m(z) G Z + as 



and 



m(z) := inf m(z, U), m(z,U) := lim max#{5„n[/|, 

U N^oon>N 



m(z) := inf m(z, U), m(z,U) := lim min =f/={S n D U\, 

— - ' u — ' — ' N 



respectively, where the infimum is taken over all open sets U containing z and S n 
is the set of poles of g n counting multiplicities. Clearly, m(z) < m(z), m(z) = oo 
if z G S, and m(z) = if z £ K. 

The forthcoming theorem implies that each pole of F attracts at least as 
many poles of meromorphic approximants as its multiplicity and not much more. 
This is one answer to question (b). 



Theorem 2.4. Let F and {<? n }ngN be as in Theorem 2.1 Then 

m(rj)>m(rj), 77 G S', (2.18) 

and 

53 (m(r ? )-m(r ? ))( 7 r-%))<F, (2.19) 

ri£S'\S 

where 

V := V(ip) + V w + (m + 2s' - 1)tt + 2 53 m ^) e ^)^ ( 2 - 20 ) 

ves'\s 

s' is the number of poles of R on S counting multiplicities, 

V w :=supV(Avg(w n ),S), (2.21) 

and 9(-) is the angle function for a system of m intervals covering S. 
We shall prove in due course that indeed Vw < +00. 

Before we proceed, we shall derive several integral representations that are 
necessary for the proofs of the above-stated theorems. To alleviate notations, it 
is convenient to formally rewrite the right-hand side of (2.3 1 as a single Cauchy 
integral. For this, we introduce for rj — x v + iy^ G C the distribution $ J; = 
\(x — x n ) ® S(y — y^), where 5 is the Dirac delta at and \ the characteristic 
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function of the non-negative semi-axis. For each k G Z + (the set of nonnegative 
integers), the partial derivative 9^ +1 <i>, ) is an analytic functional (although $ v itself 
is not) , acting on any function h holomorphic in a neighborhood of r\ according to 
the rule 

(d* +1 %,h) = hW{r,), 

where ft/ fc ) indicates the k-th derivative. Therefore, if we define A^ fc) to be d%%/k\, 
we can formally write 



dA ( v k) (t) 1 



and on rewriting R s (z) as 



Z — t (z — T]) k+1 

(n)-i 

E 

v eS' k=o 



R *( z ) = E E {z zr^k+i ■ r v,k e 



we get 



R s (z) 



dm 



where A' is given by 



m(rj) — 1 

A ' : =E E r v,k^i, k \ supp(A') = S'. (2.22) 

^GS' k=0 

This way _F can be put in the form 
with 

A:=A + A', suppA = S = S U S', (2.23) 
which makes for a convenient notation. 

Now, let {g n } be a sequence of irreducible critical points in MA(p) for some 
p G [2,oo] (cf. (pT2|). Then ([M Prop. 9.1]) 



4f (««)(£) = ln 2 £,(b n jnWn /2 ) (£) = 7« /2 (&rjn«*>n /2 ) (0- P > 2 ( 2 24) 

^(«n)(f) = I (bnU n )(0 = {b n U n y (0, P=2, 

for a.e. £ € T, where j n is some inner function, u n G -ff 2 , and h a {z) :— z^ 1 h(l/z). 
Note that if h G -ff 2 then h? G -ff 2 and vice versa. We remark that for the case 
p = 2 equation (2.24) is an interpolation condition saying that g n interpolates F 
with order 2 at the reflection of its poles. The same can be said when p G (2, oo], 
provided g n is analytic at the reflection of its poles, but this is no longer automatic 
because g n may no longer be rational. 

As usual, we denote by v n an associated singular vector to g n . According to 
(2.6), each g n can be decomposed as 
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where h n € H p . Moreover, we can write b n as q n /q n , where q n G M n and q n is the 
reciprocal polynomial of q n . Arguing like in [TOl Sec. 10] (where R s is not present), 
equation (2.241 implies easily the following orthogonality relations 



t k q n (t) 



W n (t) 



d\(t) = 0, k = 0,...,n- 1, 



(2.25) 



where w n is the outer factor of v n and A is given by (2.231. Upon rewriting (2.25) 
as 

771(77)— 1 



0. 



r]£S' k=0 

for all P n -i G V n -\ and taking P n -\ to be a multiple of Q Sl these relations yield 
for n > s 



,(t)q n (t) 



S(t) 



dA(t) = 0, k = 0,...,n 



(2.26) 



where Q s was defined in (2.4 1 



The following theorem is a result on the zero distribution of polynomials 



satisfying certain nonlinear orthogonality relations. As apparent from (2.261, it 



will be the working tool of our approach to the asymptotic behavior of poles of 
irreducible critical points, in particular of poles of best or local best approximants. 
We continue to denote by / the convex hull of S = supp(A). 

Theorem 2.5. Let {q n }n<£N be a sequence of polynomials of exact degree n with all 
zeros in D satisfying the orthogonality relations 



t k q n (t)^-d\(t) = 0, k = 0, . . . ,n - 1, 



(2.27) 



where A = A + A' is given by (2.1) and (2.22) with S C (-1,1) and S' C D. 
while W = {tL'nlri'Li * s a family of complex measurable functions on the union 
of S' and I , whose moduli are uniformly bounded above and below by positive 
constants, and whose arguments are smooth with uniformly bounded derivatives 
on I. Suppose further that A G BVT. Then the counting measures of the zeros 
of q n {z) = Ilj=i( 2 — 0,™)' namely v n :— (l/n) Ylj=i ™> converge in the weak* 
sense to ^(sj). 

The above-stated theorem is a generalization of Theorem 5.1 and Corollary 
6.2 in [5]. The main difference here is that we add a distribution of the form 



(2.221 to the measure (2.271, i.e. F may have polar singularities inside of the 



unit disk. As compared to the above reference, we simplified the proof somewhat 
by using the two-constant theorem instead of weighted potential theory. To the 
author's knowledge these are the first published results about the zero distribution 



of polynomials satisfying nonlinear orthogonality equations like (2.271 that are 



typical of rational or meromorphic approximation with free poles. In the Ph.D. 
thesis of R. Kiistner [22], an analog of [S] Theorem 5.1] is given when the measure 



Meromorphic Approximants to Cauchy Transforms 



13 



A, instead of belonging to BVT, has an argument of bounded variation and satisfies 
the so-called A-criterion introduced in [33j Sec. 4.2]: 

(L a ,• Log(l//x[* -r,t + r]) \\ fcr . 

cap teS: hmsup — w L . ' < +00 } = cap(S). 

VI r^O Log(l/r) J7 

Paralleling the arguments in [52], all the results in this section could be obtained 
under this weaker assumption, but the exposition would be heavier and we leave 
it to the interested reader to carry out the details. 



3. Proofs 

The proof of Theorem 2.5 relies on several auxiliary lemmas. 

Lemma 3.1. With the previous notation the following statements hold true 

(a) Let v be a positive measure which has infinitely many points in its support 
and assume the latter is covered by finitely many disjoint intervals: supp(y) C 
UJL-Joj, bj]. Let further ip be a function of bounded variation on supp(z^). If 
for some integer I we have 



J fl_i(t)e^ (t) di/(f) = 0, VP;_i e V1-1, 



then 

m 

52v(il>,[ aj ,bj]) > (l-m + l)7r. 
3=1 

(b) Let [a, b] C (-1, 1) and £ 6 O. Define 

g(£,t) :=Arg(i-0-2Arg(t-l/0, (3.1) 
where the term 2Arg(t — l/£) is omitted if £ = 0. Then 

V(g(Zr),la,b])<Angle(Z,[a,b}). (3.2) 

(c) Let be a real function of bounded variation on an interval [a,b], {a n (x)} a 
sequence of continuously differentiable real functions with uniformly bounded 
derivatives on [a, 6], and Q a polynomial. Then there exists a polynomial 
T^O and a constant f3 S (0, 7r/32) such that 

Arg (e t ^ {x)+a - {x]] Q{x)T{x) S j < tt/2 - 2(3 (3.3) 

for all x € [a, b] such that T{x)Q{x) ^ and all n from some infinite sequence 
NiCN. 

(d) Assume I C (0, 1) and {q n } is a sequence of polynomials of degree m n whose 
roots {£i,n, . • • , £m„,n} lie in D and satisfy 

^(tt- Angle^,/)) <C 
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where the constant C is independent of n. Then, to every e > there exists 
an integer I such that, for all n large enough, there is a polynomial T^ n of 
degree at most I satisfying: 

q n {x) 



Ti, n (x) 



< e, x e I. 



\Qn(x)\ 

In particular, the argument ofTi, n {x)/q n {x) lies in the interval (— 2e, 2e) when 
n is large enough. 

Proof, (a) This assertion follows from the proof of Lemma 3.2 in [5] if we put 
d n = I there. 

(6) When f ^ I, the proof of this statement is contained in that of Lemma 5.2 in 
[5 . In the remaining cases, one can see by inspection that (3.2 1 reduces to < ir 
when £ = b and to 7r < n when £ £ [a, b). 

(c) Observe that f(x) = ip(x) + Arg(Q(a;)) is a real function of bounded variation 
on /. Therefore by Lemma 3.4 in [5], there exist a polynomial T* =/= and a 
constant (3* € (0, 7r/16) such that 

Arg (e^^ Q(x)T* (x)\ = Arg ( e W B >T*(x)) < tt/2 - 2/3* 

for x S /, Q(x)T*(x) ^ 0. For later use we also record that, by the very construc- 
tion of T* in the cited lemma, its zeros belong to / and are discontinuity points 
of ip. 

Let K be such that |a^(x)| < K for all n € N and x € /, where the superscript 
"prime" indicates the derivative. By Jackson's theorem (cf. e.g. |28j ) there is a 
constant C > and there are polynomials {T„ ;} of degree at most I such that 

Fix / so large that CK/l < /3*/3. Being bounded of bounded degree, the sequence 
{Tn.i} has a subsequence converging uniformly on / to a polynomial Ti of degree 
at most I. Therefore, for some subsequence Ni we obtain 

/3* 



ia n (x) 



l 



Ti(x) 



< 



n e Ni, 



which implies that 



Arg (e la ^Ti{x) 



n G Ni. 



Now inequality (|3_3f follows by taking T = T*T U (3 = /3*/2, and using that 
|Arg(a +6)1 < |Arg(6)| + |Arg(6)| for any e C. 

(d) This is exactly what is proved in Lemma 5.4 of [5]. □ 



Lemma 3.2. Lei q„(z) = ]l"=i( z ~0,™) so**^ ( 2.21 ) for |^ n | < 1 /or j = 1, . . . , ra, 
w/iere A = A + A' is given by (2.1), (2.2), and (2.22), with S C (— 1, 1) and S' C D. 
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while uj n is a complex-valued measurable junction on S U S' . Consider a covering 
of S by finitely many disjoint closed intervals: S C I m := UjLi[ a ji^]- T^erJ^] 

n 

- 0&,„)) < + F(argK), 5) + £ 771(77)%) + (m + « - 1)tt, (3.4) 

where aig(u n ) is any argument function for oj n on S and m(rf) is the multiplicity 
ofrj. 

Proof. If uj n has no argument function of bounded variation on S, there is nothing 
to prove. Otherwise, we pick one and extend it to the whole of R without increasing 
the variation, as explained in Section [2] In particular, we get 

m 

F(arg(w n ), S) = V(axg(u n ), [aj,bj]). 

3=1 



As in the case of (2.261, equation (2.271 yields 



P n - s -i{t)Q s (t)q n {t)^§-d\(t) = 0, 



(3.5) 



where P„_ s _i is any polynomial in V n - s -i- 

Denote by tp n (t) an argument function for e l{p ^Q s {t)q n {t)uj n {t) /q^(t), say 

n 

ip n (t) = <p(t)+arg(«n(*))+X! m(v)^g(t-v)+J2 ( Ar §^ " &«) ~ 2Ar §^ ~ V&.n)) 

rjeS' i=l 

where it is understood that Arg(i — is omitted when ^ n = 0. It is easy to 

is of bounded ^ 

ip = ip n , dv{t) 



see that tjj n is of bounded variation. Then Lemma 3.1 'a) with 

Q s (t)q n (t)uj n (t) 



d\X\(t), and I = n 



implies that 



E V(ip n , [aj,bj]) > (n - s-m+ 1)tt. 

3=1 

So, we are left to show that 

mm n 



3 Note that the hypothesis A £ BVT is not required for this lemma to hold. 
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By the definition of ip n , we have 

m mm 

J2 V (^[a^bj]) < ^y^K^.D + ^^argK),^^,]) 



j=i 



^2^2mi V )V{Axg{'-r,),[aj,bj]) 
j=l r,eS' 
m n 

j=i »=i 



where <?(£,£) was defined in (3.1l. The assertion of the lemma now follows from 



Lemma 3.1 ^b) and the fact that, by monotonicity, V(Aig( — £), [a, b]) — Angle(£, [a, b\). 



□ 

Corollary 3.3. Let <p and arg(o; rl ) have bounded variation on S, and assume further 
that V (a,rg(u n ) , S) < C, where C is independent on n. Then, to each neighborhood 
U of S, there exists a constant kjj € N such that each q n has at most kjj zeros 
outside of U for n large enough. 

Proof. Since U is open, its intersection with (—1,1) is a countable union of in- 
tervals. By compactness, a finite number of them will cover S, say U™ =1 (aj, bj). 



Apply Lemma 3.2 to the closure of these intervals and observe that any zero of q n 



which lies outside of U will contribute to the left-hand side of (3.4 1 by more than 



some positive fixed constant which depends only on U. Since the right-hand side 



of (3.4 1 does not depend on n and is finite we can have only finitely many such 
zeros. □ 



Proof of Theorem \2.5\ Hereafter we are going to use (3.5 1 rather than (2.271 and 
we set q n (z) = Il"=i( 2 ~ 

We start by observing that we may suppose S C (0,1). Indeed, if this is 
not the case, take a negative number w such that —1 < w < a, where [a,b] = I 
denotes the convex hull of S. Then M W (S), the image of S under the Mobius 
transformation M w (z) := (z — w)/(l — zw), is a subset of (0, 1). Moreover, the 
Green equilibrium measure is invariant under Mobius transformations, i.e., for 
any Borel set E C M W (S) we have that 

^ s \E) = n%{M_ w {E)) 
(M_ w is the inverse function of M w ). This implies that the weak* convergence 
of v n to jti(5T) is equivalent to that of v™ to (j,q w ^ S \ where is the counting 
measure of the images of the zeros of q n under M w . Now, if we let 

tn{r) = q n (M- w (T))(l + WT) n , 
L s (t) = Q s (M^ w (t))(1 + wt) s , 
Pn-s-i(r) = P„- s -i(Af_ w (r))(l + U ;r)"- s - 1 , 

W n( r ) = ^n{M_ w {r)), 
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then i n is a polynomial of degree n with zeros at M t0 (^ i „), j — l,...,n. In 
addition, since M x (l/£) = l/M x (£), x £ (—1, 1), we have that 

I n (T) = q n (M- w (T))(l + WT) n . 

Analogously, L s is a polynomial of degree s with zeros at M w (rf), rj € S", and 
Pn-s-i is an arbitrary polynomial of degree at most n — s — 1. Let us show that 



satisfies orthogonality relations of type (3.5 1 for a new measure, supported this 
time in (0, 1), that still belongs to BVT. 



With the above notation equation (3.51 can be rewritten as 



s 



(P w - s - 1 Q sgra )(M-4r)) ^ ( ^-- (T \ ) e^^^)dlAl(M- w (r)) 

p„_ s _ 1 (r)L s (r)4 l (r)^ e ^*Md|A*|(r), 
m„(S) £|(r) 

where d|A*|(r) = (1+uit)g?|A|(A:/_„,(t)) is a positive measure supported on M W (S), 
<P*{t) = ip(M- w (r)) is a function of bounded variation, and {w*} is a sequence of 
measurable functions whose moduli are uniformly bounded above and below, and 
whose arguments are smooth with uniformly bounded derivatives. Further, since 
Green functions are conformally invariant, M W (S) is regular so clearly A* S BVT. 
This allows us to assume that S C (0, 1). 

First we will suppose that all zeros of the polynomials q n lie outside some 
fixed neighborhood of zero. 

For each n denote by o~ n the counting measure of zeros of q n . By the as- 
sumption that we just made there exists a compact set K such that ^ K and 
supp(V„) C K. Then supp(o-„) C K^ 1 for all n £ N, and K~ l is also compact. By 
Helly's selection theorem (cf. Thm. 0.1.3]) there exists a subsequence of natural 
numbers, Ni, such that v n — >■ v for n € Ni, where A stands for weak* convergence. 
Denote by a the reflection of v across the unit circle, i.e. da(t) = dv o (l/t). It is 
easy to check that a n A a. Observing that the assumptions on u) n imply that the 
variation of its argument on /, thus a fortiori on S, is bounded independently of 
n, it follows from Corollary |3.3| that v and a are probability measures such that 
supp(i/) C S C (0, 1) and supp(er) C S^ 1 = S^ 1 . 

Claim: it is enough to show that the logarithmic potential of v — a is a 
constant q.e. on S. Indeed, let U w ~ a — D\ q.e. on S, where 

U v ~' = J \og^—d{v-a){t). 

Then, since If 7 is harmonic outside of we have that U" is bounded quasi 
everywhere on S, hence everywhere by lower semi-continuity of potentials. Thus, 
v has finite energy and by reflection so does a. Moreover, for quasi every z G S^ 1 , 
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wc have 

U v - a {z) = 



J bg \z^i\ d[V ~ a){t) = J l0g \z-l/x\ d [{U - U) ° 

x/z 



log 



l/z 



d(a - v){x) = / log \x/z\d(a - v){x) - U v -"(l/z) 



J \og\x\d{a-u){x)-Dx=:D 2 



(3.6) 



where we used that v — a has total mass zero. Now, denote by a the balayage of 
a onto S. Then 

U 3 (t) = U a (t) + c(a;C\S) 

for quasi every t on S. Thus, as (y — <r)(C) = and since v and a have finite 
energy, we get 







£>i d(v-a)(z) 
1 



log 



\z-t\ 



U u ^{z)d{v^d){z)^ / U"-°(z)d{y-a)(z) 
d{v-a){t)d{u-a){z). 



But the energy of a signed measure is equal to zero if and only if the measure 
is zero ([32] Lemma 1.1.8]), provided that this measure is the difference of two 
positive measures with finite energy; thus, v = a. Using (3.6 1, we can obtain in 
a similar fashion that a = v, where v is the balayage of v onto S . Hence, we 
proved that a — v is the equilibrium signed measure for the condenser (S, S^ 1 ) 
( 32, Thm. VIII. 2. 6]). Then Proposition A points (b) and (e) of the appendix 
ensures that v = fi(s,T) is the Green equilibrium distribution relative to both D 
and C \ S^ 1 . Since {v n }nefii was an arbitrary weak* convergent subsequence, the 
whole sequence {v n }neN converges to Hq in the weak* sense. This proves the claim. 

Being left to prove that U V ~' J is a constant q.e. on S, suppose to the contrary 
that this is not true. Then there exist nonpolar Borel subsets of S, say E- and 
E + , and two constants d and r > such that 

U^ a (x) >d + r for x e E+ and U^^ix) <d-2r for 

In this case we claim that there exists j/o € supp(z^) such that 

U»- a (y a ) > d. (3.7) 

Indeed, otherwise we would have that 

U v [x) < U a (x) +d, xe supp(v). (3.8) 



Then the principle of domination ([321 Thm. II. 3. 2]) would yield that (3.8 1 is true 
for all z € C, but this would contradict the existence of E + . 

Since all a n are supported outside of the closed unit disk, the sequence of 
potentials {W 7 ™}^^ converges to U a uniformly on compact subsets of D. This 
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implies that for any given sequence of points {y n } C 
n — ► 00, n G Ni, we have 



lim 

+00, ngf 



U^{y n ) = U°{ m ). 



such that y„ — ► yo as 



(3.9) 



On the other hand all v n , n e Ni, have their support in D. So, by applying the 
principle of descent (j32j Thm.I.6.8]) for the above sequence {y ra }, we obtain 



liminf U v "(y n )>U v (yo). 

n — >oo, n£Ni 



Combining (3.7), ( 3.9 1 , and (3.10) we get 



liminf U^- a "(y n ) > U v 

n — >-oo n^Mi 



r M > d. 



(3.10) 



(3.11) 



Since {y n } was an arbitrary sequence in D converging to yo, we deduce from (3.11 ) 
that there exists p > such that, for any y £ [yo 
enough, the following inequality holds 



Since 



JJVn 



2p,yo 

U^-^(y)>d. 

1 Q n (y) 



2p] and n € Ni large 
(3.12) 



' (y) = - log 

n 



a n q n (y) 



where a n :— Y\.j=i \£j,n\, inequality (3.121 can be rewritten as 



Q, 



Qn(y) 



Qn(y) 



< e" 



(3.13) 



for any y € [yo — 2p, yo + 2p] and n € Ni large enough. 

Here we notice for later use that the above reasoning does not really require 
the polynomials q n to have exact degree n. Specifically, let {p n } be a sequence of 
monic polynomials of degree d n = n + o(n) where o denotes the Landau symbol 
"little oh". Moreover, suppose that the counting measures of their zeros normal- 
ized by 1 /n rather than l/d n (so it may no longer be a probability measure) are 
supported on a fixed compact set of the complex plane. Call p n these measures 
and assume that they converge to v in the weak* topology. In this case (3.10 1 and 
(3.111 still hold with v n replaced by p n , at the cost perhaps of dropping finitely 
many terms of N% and making p smaller. Thus, we obtain that 

Pn(y) 



n Qn(y) 



< e 



- nd 



(3.14) 



for any y 6 [yo — 2p, yo + 2p] and n £ Ni large enough. 

In another connection, since U u ~ a (x) < d — 2t on E. 
envelope theorem ([351 Thm. 1.6.9]) gives us 



lim inf U Vn 

n — >oo, n£Ni 



applying the lower 
(x) = U"-"[x) < d-2r, forq.e. x G E-. (3.15) 
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Let Z be a finite system of points from (—1, 1), to be specified later. Then by [2 [3] 
there existsQS'o C S such that Sq is regular, cap(i?_ n So) > and dist(Z, So) > 0, 
where dist(Z, Sq) := min z£ ^ dist(^, Sq). Put for simplicity b n (x) = q n (x) /q n (x), 
which is a finite Blaschke product. Then 

U^-^{x) = --\og\a n b n {x)\ 



and by (3.151, as \a n \ < 1, there exist N2 C Ni and x £ E_ D Sq such that 

K(x)\ > \a n b n (x)\ > e - n ^ 
for any n £ N2. Let x„ be a point where \b n \ reaches its maximum on Sq. Then 
Mn := maxima;)! = \b n (x n )\ > \a n \M n = \a n b n (x n )\ > e-^ d - T \ (3.16) 

Note that M n < 1 and therefore d — r is necessarily positive. For simplicity, we 
shall denote D \ So by D. Since the modulus of a Blaschke product is bounded by 
1 in the unit disk and log \ b n \ is a subharmonic function, the two-constant theorem 
( [3~H Thm. 4.3.7]) yields the following pointwise estimate 

log\b n {z)\<u D (z,S )log\M n \ (3.17) 

for any z £ D, where uijj(z, Sq) is the harmonic measure on D (cf. [3TJ Sec. 4.3]). 
Combining the last two inequalities we get 

K(z)\ < (M n r D ^ s °) = M n (-^j ° ( ' 0) < Mn e^ d - T ^-^^ s ^ (3.18) 

for z £ D, and obviously also when z £ Sq, where Wd(-, Sq) — 1 for Sq is regular. 
Moreover, by the regularity of Sq again, it is known ([3TJ Thm. 4.3.4]) that for any 
x £ Sq 

lim u>d{z, Sq) = 1 

z — >x 

uniformly with respect to x £ Sq. Thus, for any S > there exists r(S) < 
dist(So,T) such that for z satisfying dist(z, Sq) < r(S) we have 

1 - oj d (z, Sq) < . 



This, together with (3.181, implies that for fixed 5, to be adjusted later, we have 

K(z)\ <M n e nS , \x n -z\<r(S). 
Note that b n is analytic in D, which, in particular, yields for \z — x n \ < r(6) 

1 r b n (o 



2m J\H-x n \=r{&) (? - z ) 



7 In [5] So is n °t introduced, which makes there for a slightly incorrect argument in Theorem 3.1. 
An alternative remedy in that reference would be to apply the lower envelope theorem to Tq n 
rather than q n , as they have the same asymptotic zero distribution. 
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Thus, for any z such that \z — x n \ < r(S)/2 we get 



K(z)\ < 

Now, for any x such that 



1 4M„e™' 
2n ' r 2 (<5) 



27rr((5) 



AM n e nS 



< 



(3.19) 



the mean value theorem yields 



AM n e nS 

\b n (x) - b n (x n )\ < — \x 

r{6) 



Thus, for x satisfying (3.191, we have 



\K(x)\ > \b n (x n ) \ - \b n (x) - b n (x n ) \ > M n 



Mn 

2 ' 

Mn 

2 



2 



and, by (3.161, 



\a n b n (x) \ > 



\a n \M n 



-n(d-T) 



(3.20) 



2 - 2 

The estimates (3.141, ( 3.20[ ), together with the relation (3.51 are the main ingre- 
dients in proving the claim that XJ v ~ a is constant q.e. on S. To combine them we 
shall use a specific choice of P n _ s _i in (3.5). 



First, we pick a polynomial T such that Lemma 3.1 -(c) holds with ip = ip, 
a n = Arg(w n ), Q — Q s , and [a,b] = I for n e N3 C N2. We denote by k the 
degree of T. Second, for each n £ N3, we choose „ as in Lemma 3.1 (d) with 
e = <5/9. Note that the use of Lemma 3.1 (d) is legitimate by Lemma 3.2 and 
our assumptions on S, if and u n . Since all T/ ^ are bounded on / by definition 
and have degree at most I, which does not depend on n, there exists N 4 C N 3 
such that sequence {X^nlne^ converges uniformly to some polynomial T/ on I. 
In particular, we have that deg(Tj) < I and the argument of J](x) /q n (x) lies in 
(—(5/4,5/4) for n £ N4 large enough. Denote by 2a the smallest even integer 
strictly greater than 21 + k + s. As soon as n is large enough, since j/o G supp(^), 
there exist /3 1)71 , . . . , f3 2a , n , zeros of q n , lying in 

Of ([yo - P, Vo + p]) ■= {z e C : dist (z, [y - p, y + p]) < 7} , 

where 7, < 7 < p, will be specified later. Define 

qn&T(z)T?(z) 



The polynomial P* has degree n — s — 1 or n — s — 2, depending on the parity of 
k + s. 

Denote by /„ C (0,1) the interval defined by (3.191. By comparing (3.13) 
with (3.20), it is clear that /„ and [yo — 2p,yo + 2p\ are disjoint when n e N4 is 
large enough. 
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Now, we choose 7 = j(p) so small that 



j=l \ X Pj,nJ 



<6/2, xeR\[y -2p,y +2p]. 



J'=l 

Letting <5 be such that S < (3, the choices of T, 7] , and together imply 
Arg (p*{x)Q s {x)q n (x)'^\e i ^ 

Arg [ k„(x)| 2 • fj , ^ : • 3ftr ■ T(x)Q s (x)cj n (x)e i ^ 

< - + - + - - 2(3 < n/2 - S, 
-222 

for x € I \ [yo — 2p, 1/0 + 2p] except perhaps at points where T or Q s are equal to 
zero. This means that for such x 



Re 



\a n \\P:Q s q n ){x)^p\e^ 
Q n \X) 



> sin S 

> sin 5 



q n \x) 



a 2 MQ s TT?u n )(x) 



> 0(3.21) 



Moreover, if x G S \ [y - 2p, y + 2p] satisfies ( |3.19| , then by ( |3.16| and ( |3.20[ ) the 
above quantity is bounded from below by 

,rp ( N n / M sin< * min^g^^] |T;(x)| 2 inf ngW min xeM] |fa; n (a;)| _ 2nd+2nr 
|ilI,W,lIJI 4(diam(5) + 2p)2« 6 

= Cl |T( a; )Q s (x)| e - 2 " d + 2 -, 

where diam(S') := max Xj y 6 s |ac — y| and 

sin<5 rriin x6 [ 0i6 ] |Tj(a;)| 2 inf neN xma xe [ atb ] \u> n (x)\ 



ci := 



> 



4(diam(S') + 2p) 2a 

by the construction of 7] and the uniform boundedness of {|w n |} from below. Thus 



Re 



> sin S 



\a n \ 2 P:(t)Q s (t) qn (t)^e^dm) 

S\[y -2p,y +2p\ InK 1 ) . 



S\[y a -2p,y +2p] 



alP:(t)Q s (t) qn (t)^-e^ 

<ln\ l ) 



d\X\(t) 



> c ie - 2nd+2nT \T(t)Q s (t)\d\X\(t)>c 2 e- 2nd+ni2T - L5 l (3.22) 

Jsm n 

The last inequality is true by the following argument. First observe that from 
( 3.20 1 and 3.13 ) that J„ n [yo — 2p, y + 2p] = for all n large enough. Next, recall 
that x n , the middle point of /„, belongs to Sq, where dist(So, Z) > and Z is 
a finite system of points that we now choose to be the zeros of TQ S in (—1, 1) if 
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any. Then TQ S , which is independent of n, is uniformly bounded below on /„ and 



(3.22 



and 



follo ws from this and hypothesis BVT point (2). On the other hand (3.131, 
3.141 applied with p n = P*/ (leading coefficient of TTf), yield that 



[i/o-2p,yo+2p] 



d\X\(t) 



< c 3 e" 



-2nd 



(3.23) 



where we used uniform boundedness of {|w n |} from above. Choosing 5 so small 
that 2r — L5 > 0, which is possible, the bound in (3.221 becomes bigger than the 



bound in (3.23) for n large enough. But this is impossible, since by (3.5) the sum 



of these two integrals must be zero. 

We just completed the case when all the zeros of polynomials q n stay away 
from the point zero. Now we shall consider the general situation. Let e > be 
such that U :— D \ {z : \z\ < e} is a neighborhood of S. Corollary 3.3 says that 



there exists a constant kjj G N such that each q n has at most kjj zeros outside of 
U , that is zeros which have modulus less than e. In this case, from any sequence 
of natural numbers, we can extract a subsequence, say No, such that for some 



number m < kjj, 
these zeros £i n , . 



q n has exactly m zeros outside of U for each n G No- Denote 
• )£mn; an d consider the polynomials 



q n (z) 



n G N . 



n£i(i-^,j' 

Then the sequence {q n }n<£N will satisfy the following weighted orthogonality re- 
lation: 



Pn- 



.i(t)Q s (t)q n (t) 7 ^r,dX(t) = 0, 



feW) 5 



where 



Pn- 



tel. 



G V n 



In what follows we are going to stress the modifications needed to adapt the 
previous proof to the present case. Let, as before, Ni C No be a subsequence of 
natural numbers such that v n — * v, n G Ni. Because we only discarded a fixed 
number of zeros from q n (z) to obtain q n , the counting measures cr* of the zeros 
of q* n (normalized by 1/n), again converge weak* to a. Since If 7 " enjoys all the 
relevant properties of ?7 <T ™, inequalities ( |3.13 1 and (3.141 remain valid with q n 
replaced by q n . 

Further, define b n (x) as q n {x)/q n {x). In this case, b n is no longer a Blaschke 
product, but rather a Blaschke product multiplied by the polynomials™^ (z— £j, n )- 



Then we get instead of (3.17) that 



log|M*)l < u D (z,S )log\M n \ + (1 - u D (z, S ))mlog2, 



and (3.181 can be replaced by 



|6 n (z)| < M ne n ( d - T + mlo e 2 / n ){ 1 ~^n(z,S )) < ^ ne n{d-T+l)(l-w D {z,S Q )) ^ 
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for n large enough. This yields an insignificant modification of r{$) (we should 
make 1 — u>d{z, So) less than S/(d — t + 1) rather than just 5 /{d — t)). Lemma 
|3.1[ d) can be applied to the polynomials g* rather than q n without change. 

Therefore we are left to show that {o>* } is uniformly bounded above and below 
on /, and that its arguments are smooth with uniformly bounded derivatives on 
I with respect to n. The uniform boundedness of {w*} easily follows from the 
estimates 



2m 



< 



n 



i 



(i - **,,„) 



< 



1 - e 



2 m 



Since none of the £j >n , j = 1, . . . ,m can come close to / 1 , Arg(l — t£j n ) is 
a smooth function on / whose derivative Im^, n /(1 — i£ ■ n ) is uniformly bounded 
there independently of j and n. Then the rest of the assumptions on {w*} follows 
from the representation 



arg«(i)) = arg(w„(t)) 
This completes the proof of the theorem. 



2^Arg(l-^- 

.7=1 



□ 



The forthcoming lemma is needed for the proof of Theorem 2.1 Recall that 
H, a family of functions analytic in some fixed domain of the complex plane, is 
called normal if each sequence of functions from TL contains a locally uniformly 
convergent subsequence. 

Lemma 3.4. Let p € (2, oo] and {g n }neN be a sequence of a sequence of irreducible 
critical points of order n o/MA(p) for F given by (2.3) and (2.4-) with A sa£zs/?/m(j^] 
(2.1) and (2.2). Further, let v n be an associated singular vector to g n with inner- 
outer factorization given by v n — b n ■ w n for some Blaschke product b n and w n 
an outer function in H p . Then the families W := {w n } and W p < := are 

normal in D and C \ S* respectively, where S* denotes the reflection of S across 
T. Moreover, any limit point of ' W is zero free in D. 

Proof. The main idea of the proof was given in [101 Thm. 10.1]. The necessary 
modification for the case of complex measures with argument of bounded variation 
were given in [7J Prop. 6.3]. Only a simple adjustment is needed in the present case 
where the approximated function may have polar singularities inside D. Namely, 
after the initial choice of a polynomial T in j7j Prop. 6.3, Eq. 6.12] has been made, 
one systematically applies the arguments in this proposition with TQ 2 instead of 
T. Another modification that one has to introduce is to define G n in 10, Thm. 
10.1] as 

1/2 

G n (z) := (b n Qw n ' )(z), 



8 Notc that we do not require the hypothesis A £ BVT to hold. It is sufficient for the lemma to 
hold to have a measure with an argument of bounded variation and infinitely many points in the 
support. 
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in order to obtain the desired bound in [7j Eq. 6.13]. The interested reader can 
find the full proof of this lemma in 34, Lem. 3.11]. □ 

Note that in the previous lemma the normality of W in D was clear before- 
hand by the Cauchy formula, since ||iy n ||p' — 1, but the nonzeroing of every limit 
point was not. 



Proof of Theorem \2.1\ Let v n be a singular vector associated to g n with inner- 
outer factorization given by v n — b n ■ w n for some Blaschke product b n — q n /q n 
and some outer function w n € H p , \\w n \\ p i = 1, where w n = 1 when p = 2. The 
poles of g n are exactly the zeros of q n . Moreover, {q n } is a sequence of polyno- 



mials satisfying weighted non-Hcrmitian orthogonality relations (2.25). Thus, the 
assertion of the theorem will follow from Theorem |2.5| if only VV = {w n } is uni- 
formly bounded above and below on /, the convex hull of S, and if it is a family of 
functions whose arguments are smooth with uniformly bounded derivatives on /. 
In the case p — 2 this is trivial since each w n = 1. In the case p € (2, oo] Lemma 



3.4 says that VV is a normal family. Thus, it is uniformly bounded above on /. 
Moreover, since all limit points of VV are zero free in D, this family is uniformly 
bounded below on I (in fact on any compact subset of D). Further, the derivatives 
again form a normal family and so does the logarithmic derivative w' n /w n in D. 
Since the imaginary part of the latter is equal to d a,rg(w n ) / dt on /, we see that 
the desired conditions on VV are satisfied. □ 

Before we prove Theorem |2.2| we shall need one auxiliary lemma. 

Lemma 3.5. Let D be a domain in C with non-polar boundary, K' be a compact 
set in D, and {u n } be a sequence of subharmonic functions in D such that 

U n {z) < M - e n , z e D, 

for some constant M and a sequence {e n } of positive numbers decaying to zero. 
Further, assume that there exist a compact set K' and positive constants e' and 5' , 
independent of n, for which holds 

u n (z) < M — e', z G K n C K\ cap(K n ) > 5'. 

Then for any compact set K C D \ K' there exists a positive constant e(K) such 
that 

u n {z) <M — e(K), zeK, 

for all n large enough. 

Proof. Let uo n be the harmonic measure for D n := D\K„. Then the two-constant 
theorem (3T] Thm. 4.3.7] yields that 

u n {z) < {M-e')uj n (z,K n ) + (M-e n )(l-u n (z,K n )) 
< M - (e' - e n )iv n (z,K n ), z e D n . 

Thus, we need to show that for any K C D\K' there exists a constant S(K) > 
such that 

u n (z,K n )>8(K), zeK. 
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Assume to the contrary that there exists a sequence of points {z n }neNi C K, 
Nj C N, such that 



co n (z n , K n ) — > as n — > oo, n E N%. 



(3.24) 



By [3TJ Theorem 4.3.4], u> n (-, K n ) is the unique bounded harmonic function in D n 
such that 

lim u;(z,K n ) = l/f„(C) 

for any regular ( E dD n , where ljc n is the characteristic function of dK n . Then it 
follows from (A. 8 1 of the appendix that 

cMK n ,dD)U^ 0D) 



w n (',K n ), 



(3.25) 



where (J-(K n ,dD) is the Green equilibrium measure on K n relative to D. Since all the 
measures H(K n ,dD) are supported in the compact set K', there exists a probability 
measure \i such that 

^(K n ,dD) — > A 1 a s n ~> co, n £ N2 C Ni. 

Without loss of generality we may suppose that z n — > z* € if as n ~> 00, n € N2. 
Let, as usual, go{-,t) be the Green function for D with pole at t E D. Then, by 
the uniform equicontinuity of {<?£>(•, t)}t^K' on if, we get 



U£(z*)^0 as n->oo, n e N 2 . 



Therefore, (3.241 and ( 3.25| ) necessarily mean that 

cap(if„, dD) ^0 as n — > cxd, 71 e N2. 



(3.26) 



By definition, l/cap(if„, dD) is the minimum among Green energies of probability 
measures supported on K n . Thus, the sequence of Green energies of the logarithmic 
equilibrium measures on K n , (J-K n , diverges to infinity by (3.261. Moreover, since 

{s(-,t)+lGg|.-t|} teJr , 

is a family of harmonic functions in D whose moduli are uniformly bounded above 
on K' , the logarithmic energies of [ijc n diverge to infinity. In other words, 

cap(if n ) ^0 as n— > 00, n G N2, 

which is impossible by the initial assumptions. This proves the lemma. □ 



Proof of Theorem 2.2 To prove the convergence in capacity we first establish an 
integral representation for the error [F — g n ). As usual we denote by v n = b n w n 
a singular vector associated to g n , where b n is a Blaschke product of degree n 
and w n is an outer function. By (2.111 and Fubini-Tonelli's theorem, we have for 
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z e B\S 

r n (z) 



(F-g n )(z) = 
1 



A F {v n ){z) = P_{Fv n )(z) 
v n (z) v n {z) 



1 

v n (z) Jt z-£ 27TZ v n (z) 

1 fvn{t) d\(t) = Uz) 



( Z -m-t) dX{t) ^i 



q n (t) w n (t) 

v n (z) J z-t~"' y "' q n (z)w n (z) J z-t q n (t) 
In another connection, the orthogonality relations ( |2.25 1 yield 

J z-t qi(t) 
This, in turn, implies that ( |3.27 1 can be rewritten as 

qt( z ) f Qn(t) W n (t) 



dX(t). (3.27) 



r n {z) 



q n (z)w n {z) J z-t q%(t) 



d\(t), zeO\S. 



(3.28) 



Since the majority of the zeros of q n approach S by Corollary |3.3| we always can 
choose s of them, say £i,n, ■ ■ ■ , £s,n, in such a manner that the absolute values of 



h,n{z) ■= - £j>) and l s>n (z) := z s l s , n (l/z) 



(3.29) 



are uniformly bounded above and below on compact subsets of C \ S and in D, 
respectively, for all n large enough. Using orthogonality relations (2.251 once again 
and since Q s vanishes at each r\ € S' with multiplicity m(?y), we can rewrite (3.28) 
as 

(&(z) [ (q n q* n Q s w n )(t) d\(t) 



r n (z) 



(q n qnQsW n )(z) 



z-t 



where q*(z) := q n (z) /l s , n (z). Set 

(\q* n \ 2 QJs,nW n )(t) d\(t) 
cfc(t) z-t' 

(qnqnQsWn)(z)' 



B n (z) := J 



zeC\S, 



(3.30) 



(3.31) 



so that 



r n {z) 



First, we show that 



| Bn |i/2«- p exp{-l/cap(S,T)} 



(3.32) 



on compact subsets of C \ S. Denote, as usual, b n — q n /q n . Then for any compact 
set K c C \ S there exists a constant c(K), independent of n, such that 



\B n {z)\ <c{K)\\b n f s , zeK, 



(3.33) 
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by the choice of l Sjn and Lemma 3.4 Let v n be the counting measures of zeros of 
b n . Then 

limsup |6 n (t)| 1 ^" = limsupexpj— U^"(t)} — exp{— liminf U^"(t)} 

n — >oo 

= exp {-Uo (s ' T} (t)} = exp{-l/cap(S*,T))} for q.e. t(g.m) 



by Theorem 2.1 the lower envelope theorem (351 Thm. 1.6.9], and (A. 8 1 of the 
appendix. Moreover, by the principle of descent (3H Thm. 1.6.8], we get that 

limsup|6„(i)| 1/Tl < exp{-l/cap(5,T)} (3.35) 

n — >oo 

uniformly on S. It is immediate from (3.341 and ( |3.35 I that, in fact, 



lim 116, 



"US 



exp{-l/cap(S*,T)}. 



(3.36) 



Suppose now that (3.32 1 is false. Then there would exist a compact set K' C C \ S 
and e' > such that 



capjzeif': \B n {z)\ 1 ' 2n - exp{-l/cap(S*, T)} > e'} ■/* 



0. 



(3.37) 



Combining (3.371, (3.36), and ( |3.33 ) we see that there would exist a sequence of 
compact sets K n C K' , cap(K n ) > 6' > 0, such that 



\B n {z)\ l ' 2n < exp{-l/cap(5,T)} - e', z e K„ 



(3.38) 



for all n large enough. Now, let L be a closed Jordan curve that separates S from 
K' and contains S in the bounded component of its complement. Observe that 
(l/2n) log \B n \ is a subharmonic function in C \ S. Then ( |3.33 1, (3.361, and (3.38) 
enable us to apply Lemma 3.5 with M = — l/cap(S f , T) which yields that there 
exists e(L) > such that 



\B n (z)\^ 2n < exp{-l/cap(S,T) - e(L)} 
uniformly on L and for all n large enough. Define 

dz 



(3.39) 



J 71 ■ — 



Tf(z)T(z)l s , n {z)B n {z) 



2iri 



where the polynomials T} and T are chosen as in Theorem 2.1 (see discussion 
after ( 3.20 >) . Then if the limit in (3.32) would not hold, we would get (3.391 and 
subsequently 



limsup Jl /2n < exp{-l/cap(S*,T) - e(L)}. 
In another connection, the Cauchy integral formula yields that 

7 f T*( W ^ 7— FT ( f (K\ 2 Qsls,nW n )(t) dX(t) \ dz 

./„ = I T t (z)T(z)l s , n (z) [J m — t )^- 

\ql{t)\3^{TQ s w n ){t)e^d\\\{t) 



(3.40) 



(3.41) 
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Exactly as in (3.21), we can write 
Re 



(JfTQ„w n )(t)e^ 



> sin(<5) 



(T t TQ s w n ){t) 



t e i, 



(3.42) 



where / is the convex hull of S and <5 > has the same meaning as in Theorem 



2.1 (see construction after (3.18)). Thus, we derive from (3.411 and (3.421 that 

(3.43) 



J n >sm(S) J \b 2 n (t)\ \{T?TQ s w n ){t)\d\\\{t). 

Let So be a closed subset of S of positive capacity that lies at positive distance 
from the zeros of TQ S on / (see Theorem 2.1 for the existence of this set). Further, 
let x n € S be such that 

HMSo = \bn(x n )\- 

It follows from (3.341 and (3.361 that 



lira 1 1 b n 



I So 



exp{— l/cap(5, T)}, 



and therefore 



\\b n \\ So > exp{-n(e+ l/cap(5,T))} 



for any e > and all n large enough. Proceeding as in Theorem 2.1 (see equations 

(3.44) 



(|3.19| and (|3.20|), we get that 



|&n(*)l > -exp{-n(e + l/cap(S,T))}, t e 7„, 



where 



I n ■= {x G S : |i — x n | < r^e n5 } 
and r$ is some function of 5 continuous and vanishing at zero. Then by combining 
(3.431 and (3.441, we obtain exactly as in (3.221 that there exists a constant C\ 
independent of n such that 



Jn > sin(<5) J \b n {t)\ \(TfTQ s w n ){t)\d\X\(t) > c t exp{-2n(e+£5/2+l/cap(S,T))}. 

Thus, we have that 

liminf jy 2n >exp{-e-L6/2- l/cap(5,T)}. (3.45) 

n — *oo 

Now, by choosing e and S so small that e + LS < e(T), we arrive at contradiction 
between (3.40) and (3.451. Therefore, the convergence in (3.321 holds. 

l/2n 



Second, we show that 



cap 



exp{U^(z)} 



(3.46) 



q n (z)q n {z)Qs(z)w n (z) 

on compact subsets of D \ S. Let K C D \ S be compact and let U be a bounded 
conjugate-symmetric open set containing K and not intersecting S. Define 



bn,i(^) := 



n 



Cet/: 6„(C)=0 



and 6„, 2 (z) := b n (z)/b n . 1 (z). 
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Further, let q n ,\ and q n .2 be the numerators of b n ,\ and 6 n ,2, respectively. Corollary 

ch 6„,i : 



373] yields that there exists fixed m £ N such that each b n ,i has at most m zeros. 

l/n 



Then 



1 



|& n , 2 W| 1/n ->exp{-< (s ' T) (z)} and 

Qn,2\Z) 

uniformly on K by Theorem |2.1| Moreover, it is an immediate consequence of the 
choice of l s , n , the normality of {w n } in D (Lemma 3.4 1, the uniform boundedness 
of the number of zeros of Qn,i(z)q nt i(z)Q a (z), and the lemniscate theorem [3"T1 
Thm. 5.2.5] that 

l/2n 



qiAz)l Stn (z) 



qnAz)q n ,i{z)Q s {z)w n (z) 



cap 



1, z£ K. 



Thus, we obtain (3.46). It is clear now that (2.13) follows from (3.301, (3.321, and 

Let us finally fix p = 2 and p = oo . In this former case w n 
and in the latter {w n } is a normal family in C \ S* by Lemma 



3.4 



1 for any n £ N 
and therefore 



r n (z) is de fi ned ev e rywhe re out side o f S U S*. The limit in ( |2.14 1 easily follows 
from ( |2.13[ ), ( |3.30|) ( |3.32| , and ( |3.46[ ) since |6 n (z)| = |6„(l/z)pln_C. It remains 
only to show (2.151 for p = 2 (in this case the error is defined in C \ S). Taking 
into account (3.361, (3.33), and the above-mentioned symmetry, it is sufficient to 
prove that 

limsup|6 n (z)| 1 /"<exp{-f/; (s ' T) (z)} 



uniformly in D. The latter follows by Theorem |2.1| the principle of descent, and the 
following fact. Let a £ (0, 1) be such that ScD a and denote by b n , a the Blaschkc 
product with those zeros of b n that are contained in D a . Then |6„| < \b n a \ in D 
and the zeros of b n ^ a have the same limiting distribution (/^(s,t)) as the zeros of b n 
by Corollary 3.3 Thus, Ifon.al 1 ^™ converge locally uniformly in {z : a < \z\ < l/a} 
to some function that coincides with exp {-< <S ' T) (^)} on {z : a < \z\ < 1}. This 
finishes the proof of the theorem. □ 



Proof of Corollary \2.3\ Let {r n } be a sequence of solutions of MA(2) (see (2.7l). 
Then {r n } is a sequence of rational functions with poles in D for which (2.15) 
holds. Thus, 



limsupp„(F,E) 1/2n < limsupHF- 



|l/2r, 



' n Hi 



< 



exp 



1 

"cap(5,T) 



since F — r n are analytic on E. On the other hand, AAK theory implies that 

fflloo] =liminfcr n (A F ) 1 / 2n . 



lim Mp n (F,E) 1/2n > liminf inf ||F 

n — >oo n — ^oo Vp£-^^° 



It follows from (2.11 1 and (2.241 that for every g n , best AAK approximant of order 
n, it holds that 

\F - g n \ = o- n (A F ) a.e. on T. 
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Therefore, we get from (3.30) that 

<7n(0 ls,n(0Bn(0 



<Jn{A F ) 



<?«(0 Q.(£K(0 



for a.e. (eT, 



where L 



and B n were defined in (3.29) and (3.311, respectively. Since 

l/r 



on compact subsets of 1 
and Lemma [3^41 that 



lgn(g)j 
l<7n(z)| 

\ S by Theorem 



cap 



2.1 



we immediately deduce from (3.32) 



lim <J n (A F ) 



l/2n _ 



exp 



1 

"cap(5, T) 



which finishes the proof of the corollary. 



□ 



Proof of Theorem\2.4\ Inequality (2.181 is trivial for any 77 e S'nS. Suppose now 



that 1] S S' \ S and that 777(77) < 777(77). This would mean that there exists an 
open set U, U n S — {if}, such that 777(77, U) < 771(77) an d therefore would exist a 
subsequence Ni C N such that 

#{S n nU} <m{r]), ueNi. 



It was proved in Theorem 2.2 that {g n } converges in capacity on compact subsets 
of D\ S to F. Thus, {<7n}raeNi is a sequence of meromorphic functions in U with at 
most 777,(77) poles there, which converges in capacity on U to a meromorphic function 
F\jj with exactly one pole of multiplicity 777(77). Then by Gonchar's lemma |16[ 
Lemma 1] each g n has exactly 771(77) poles in U and these poles converge to 77. This 
finishes the proof of (2.181. 



3.3 



Now, for any r\ G S' \ S the upper characteristic 777(77) is finite by Corollary 
Therefore there exist domains D v , D v n S = {77}, such that 777(77) = ™>{v> D^), 
77 <E S' \ S. Further, let O(-) be the angle function defined in (2.17) for a system of 
to intervals covering S and let S n — ■ • ■ , £n,n}- Then by Lemma 3.2 we have 



U - 0&>)) < V(ip) + V w + (m + a - 1)tt + ^ ^(0^(0, 



(3.47) 



where Fyy was defined in (2.211. The finitencss of Vyy was obtained in the proof 



of Theorem 2.1 Then for 77 large enough (3.47) yields 



E E fr-e&,*))-Mv)(*-8(ri))]<v, 

ves'\s V^e-D,, 
where V was defined in ( |2.20 1. Thus, 



E (#{S n n D v }- 777(77)) (tt-%)) < E #{SnnI^}Uax0(O-%) 



+v 



(3.48) 
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for all n large enough. However, since {max n >jv #{5* n n -D^jjjveN is a decreasing 
sequence of integers, 777(77) = fn(rj,D v ) — #{S n n D v } for infinitely many n E N. 



Therefore, we get from (3.48) that 



neS'\s 



(771(77) — 771(77)) (it — 0(rf)) < V - 



neS'\s 



771(77) ( max 0(£) 



6(7))). (3.49) 



Observe now that the left-hand side and the first summand on the right-hand side 



of (3.491 are simply constants. Moreover, the second summand on the right-hand 



side of (3.491 can be made arbitrarily small by taking smaller neighborhoods D v . 
Thus, ([2TT9F follows. □ 



4. Numerical Experiments 

The Hankel operator Ap with symbol F £ H°° + C(T) is of finite rank if and 
only if F is a rational function |241 Thm. 3.11]. In practice one can only compute 
with finite rank operators, due to the necessity of ordering the singular values, so a 
preliminary rational approximation to F is needed when the latter is not rational. 
One way to handle this problem is to truncate the Fourier series of F at some 
high order N. This provides us with a rational function Fn that approximates 
F in the Wiener norm which, in particular, dominates any LP norm on the unit 
circle, p G [l,oo]. It was proved in |21j that the best approximation operator 



from (mapping F to g n according to (2.10 1) is continuous in the Wiener 
norm provided (n + l)-st singular value of the Hankel operator is simple. It was 
shown in [6, Cor. 2] that the last assertion is satisfied for Hankel operators with 
symbols in some open dense subset of H°° + C(T), and the same technique can 



be used to prove that it is also the case for the particular subclass (2.3 1. Thus, 
even though the simplicity of singular values cannot be asserted beforehand, it is 
generically true. When it prevails, one can approximates Fn instead of F and get 
a close approximation to g n when N is large enough. This amounts to perform the 
singular value decomposition of Ap N (see [33 Ch. 16]). When 2 < p < 00 there 
is no difficulty with continuity issues, but the computation of g n has to rely on a 
numerical search. To numerically construct rational approximants when p = 2, we 
used the above truncation technique together with the Hyperion software described 
in ED]. 

In the numerical experiments below we approximate function F given by the 
formula 

J[-6/7,-l/8] z ~ t J [2/5,1/2] t — 2i z — t J [2/3,7/6] z ^ t 

2 6 24 



(z + 3/7-4i/7) 2 (z-5/9-3i/4) 3 (z + 1/5 + 6i/7) 4 ' 
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FIGURE 1. AAK (left) and rational (right) approximants to F of de- 
gree 8 

On the figures the solid lines stand for the support of the measure, diamonds 
depict the polar singularities of F, and circles denote the poles of the corresponding 
approximants. Note that the poles of F seem to attract the singularities first. 



Appendix 

Below we give a brief account of logarithmic potential theory that was used exten- 
sively throughout the paper. We refer the reader to the monographs [311 EZj for a 
complete treatment. 



34 



L.Baratchart and M.Yattsclcv 




0.8-0.6-0.4-0.2 0.2 0.4 0.6 0. 



FIGURE 2. AAK (left) and rational (right) approximants to F of de- 
gree 13 



The logarithmic potential and the logarithmic energy of a finite positive mea- 
sure /i, compactly supported in C, are defined by 



U"{z) := J log ^—dnit), zeC, 



(A.l) 



and 



IM ■■= J U"(z)d»(z) = J J togj-l^dM*)dM*), (A.2) 

respectively. The function Z7 M is superharmonic with values in (—00, +00], which 
is not identically +00. It is bounded below on supp(^i) so that J[/z] G (—00, +00]. 
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Let now E C C be compact and A(E) denote the set of all probability 
measures supported on E. If the logarithmic energy of every measure in A(E) is 
infinite, we say that E is polar. Otherwise, there exists a unique he G A(E) that 
minimizes the logarithmic energy over all measures in A(E). This measure is called 
the equilibrium distribution on E. The logarithmic capacity, or simply the capacity, 
of E is defined as 

cap(-E) = exp{-I[ii E ]}. 
By definition, the capacity of an arbitrary subset of C is the supremum of the 
capacities of its compact subsets. We agree that the capacity of a polar set is zero. 
We say that a sequence of functions {h n } converges in capacity to a function h on 
a compact set K if for any e > it holds that 

lim cap({z £ K : \h n (z) - h(z)\ > e}) = 0. 

n — >oo 

Another important concept is the regularity of a compact set. We restrict to 
the case when E has connected complement. Let g^ E (-,t) be the Green function 

of C \ E with pole at t £ C \ E, i.e. the unique function such that 

(i) g^ E (z, t) is a positive harmonic function in (C \ E) \ {t}, which is bounded 
outside each neighborhood of t; 

(ii) 9c\e( z > *) ~ { _ bg |' z - t | } J ^ is bounded near *! 

(iii) lim gpx E (z,t) = for quasi every £ € E. 

Points of continuity of 3c\b(''^) on ^eE, the outer boundary of E, are called 
regular, other points on d e E are called irregular; the latter form a polar set. If 
every point of d e E is regular, we say that the whole set E is regular. 

Throughout we use the concept of balayage of a measure (|32j Sec. II.4]). Let 
I) be a domain (connected open set) with compact boundary dD whose comple- 
ment has positive capacity, and /i be a finite Borel measure with compact support 
in D. Then there exists a unique Borel measure fi supported on dD, with total 
mass is equal to that of /j,: ||/i|| = ||/i||, whose potential i7 M is bounded on dD and 
satisfies for some constant c(/i; D) 

UV(z) = £/ M (z) +c(/j,;D) for q.e. zeC\D. (A.3) 

Necessarily then, we have that c(/x; D) — if D is bounded and c(/i; D) = 



J gn(t,oo)dfi(t) otherwise. Equality in (A.3 1 holds for all z G C \ D and also 



at all regular points of dD. The measure fi is called the balayage of fi onto dD. It 
has the property that 

U^(z) < U»(z) + c(fi; D) for every z e C, (A.4) 

and also that 

hdjjL= I hdjl (A.5) 



for any function h which is harmonic in D and continuous in D (including at infinity 
if D is unbounded). From its defining properties ~j2 has finite energy, therefore 
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it cannot charge polar sets. Consequently, on solving the generalized Dirichlct 
problem |31l Thm. 4.1.5] for an arbitrary positive continuous function on 3D, 



it follows from (A. 5) that the balayage of a probability measure is a probability 
measure. 

The minimal energy problem can also be formulated for signed measures 
[3"2l Thm. VIII. 1.4]. In particular for E\, E 2 two disjoint compact sets of positive 
capacity, there exists a unique measure fi* = fj,* — /i^, with fj,* E A(£i) and 
/i| G A(E 2 ), that minimizes the energy integral 

- M2] = J log j^—f^i ~ ^2){t)d(m - // 2 )(z), fij G A(^), j = 1,2. 

(A.6) 

It can be proved ([32l Lemma 1.8]) that I\fa"} is positive and finite. The value 
cap(£'i, E 2 ) = l/I[/J,*} is called the condenser capacity of the pair (E 1 ,E 2 ). Further, 
it holds that [i\ = fa and ^2 = Mi> where fa (resp. fa) indicates the balayage of 
^1 (resp. fi 2 ) onto d(C \ E 2 ) (resp. d(C \ Ei)); this property in fact characterizes 
fi*, see [311 Thm. VIII.2.6]. 

In analogy to the logarithmic case, one can define the Green potential and the 
Green energy of a positive measure /i supported in a domain D with compact non- 



polar boundary. The only difference is now that, in (A.1|-(A.2), the logarithmic 
kernel log(l/|z — t\) gets replaced by gi)(z,t), the Green function for D with pole 
at t G D. The Green potential relative to the domain D of a finite positive measure 
/i compactly supported in D is given by 



U£(z) = j g D (z,t)dfi(t). 

It can be re-expressed in terms of the logarithmic potentials of [i and of its balayage 
ju onto dD by the formula [32, Thm. II.4.7 and Thm. II. 5.1] 

U^(z)=c{ t i;D)-U^{z), zeD, (A.7) 



where c(/Lt; D) was defined after equation (A.3l. Moreover, (A.7 1 continues to hold 
at every regular point of 3D; in particular, it holds q.e. on dD. 

Exactly as in the logarithmic case, if E is a compact nonpolar subset of D, 
there exists a unique measure H(e,8D) S A(_E) that minimizes the Green energy 
among all measures in A(E). This measure is called the Green equilibrium distri- 



bution on E relative to D. By (A.7 1 we have that 

Up {E ' 8D) (z) = U^ E - 3 ^{z)-U^^ r ^{z)+c{^ (E , dD) ;D), ZED, and q.e. z G 3D, 

where faj^QD) is the balayage of ^(e,3D) onto dD. In addition, the Green equilib- 
rium distribution satisfies 

^(,) = _^ ) for q.e. ZEE, (A.8) 

where cap(-E, dD) is Green (condenser) capacity of E relative to D which is the 
reciprocal of the minimal Green energy among all measures in A(E). Moreover, 



equality in (A.8 1 holds at all regular points of E. 
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For the reader's convenience, we formulate below a proposition that was of 
particular use to us. It has to do with the specific geometry of the disk, and we 
could not find an appropriate reference for it in the literature. The proof of this 
proposition can be found in [34, Prop. A.l]. 

Proposition A Let E C D be a compact set of positive capacity not containing 
with connected complement, and E* stand for its reflection across the unit cir- 
cle, i.e. E* := {z £ C : 1/z £ E}. Further, let ji £ A(E) and a £ A(E*) solve the 
signed energy problem for the condenser (E,E*). Then, we have that 

(a) a is reflected from fx across the unit circle, i.e. <j{B) — fx(B*) for any Borel 
set B, and likewise \x is reflected from er; 

(b) \x is the Green equilibrium distribution on E relative to C \ E* and a is the 
Green equilibrium distribution on E* relative to C \ E; 

(c) fx — a, where A denotes the balayage of the measure A on T. Moreover, the 
balayage of Jx onto E is fx and the balayage of fx onto E* is a; 

(d) fx is the Green equilibrium distribution on T relative to both C\E and C \ E* ; 

(e) fx is the Green equilibrium distribution on E relative to D and a is the Green 
equilibrium distribution on E* relative to C \D). 
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